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CALCULATIONS  OF  OPTIMUM  OPERATING  PARAMETERS 
FOR  A GYROTRON  TRAVELLING  WAVE  AMPLIFIER* 


I.  INTRODUCTION 

The  gyrotron  is  a new  type  of  microwave  device  employing  the  electron  cyclotron  maser 
mechanism.  It  ideally  consists  of  an  ensemble  of  monoenergetic  electrons  following  helical 
trtjectories  around  the  lines  of  an  axial  magnetic  field  inside  a fast  wave  structure  such  as  a 
metallic  tube  or  waveguide.  The  physical  mechanism  responsible  for  the  radiation  in  the  gyro- 
trons  has  its  origin  in  a relativistic  effect.  Initially,  the  phases  of  the  electrons  in  their  cyclo- 
tron orbits  are  random,  but  phase  bunching  can  occur  because  of  the  dependence  of  electror 
cyclotron  frequency  on  the  relativistic  electron  mass.  Those  electrons  that  lose  energy  to  the 
wave  become  lighter,  rotate  faster,  and  hence  accumulate  phase  lead  while  those  electrons  that 
gain  energy  from  the  wave  become  heavier,  rotate  slower,  and  accumulate  phase  lag.  This  can 
result  in  phase  bunching  such  that  the  electrons  radiate  coherently  and  amplify  the  wave.  En- 
ergy transfer  from  the  electrons  to  the  wave  is  optimized  when  <u  —k.  v.  — 0,  where 

(o,k^  v^  s,  {If,  are  respectively,  the  wave  frequency,  axial  wave  number,  axial  electron  velocity, 
cyclotron  harmonic  number  and  relativistic  electron  cyclotron  frequency.  Early  descriptions  of 
the  physical  process  are  to  be  found  in  the  works  of  Twiss ',  Schneider^,  and  Gapanov^,  while 
the  first  deliberate  experimental  study  was  made  in  1964  by  Hirshfield  and  Wachtel^.  An  ex- 
tensive description  of  the  early  work  is  to  be  found  in  recent  review  papers*'®. 


•Work  supported  in  part  by  the  Naval  Electronic  Systems  Command.  Task  XF  54-587,  by  the  Army  Ballistic  Missile 
Defense  Advanced  Technology  Center,  MIPR  W31RPD-73-Z787,  and  by  the  Naval  Surface  Weapons  Center 
(Dahlgren),  Task  SF32-302-4IB. 
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The  gyrotron  emits  radiation  near  the  frequency  <u  - fl*  + so  the  wavelength  is 
primarily  determined  by  the  strength  of  the  applied  magnetic  held  and  not  by  the  dimensions 
of  a resonant  structure.  Thus  unlike  most  other  microwave  tubes,  the  internal  dimensions  of 
the  device  may  possibly  be  large  compared  to  the  wavelength,  and  high  power  handling  capa- 
bility becomes  compatible  with  operation  at  millimeter  and  submiliimeter  wavelengths. 
Indeed,  the  highest  recorded  millimeter  wave  power,  both  peak  and  average,  have  been  at- 
tained in  gyrotron  devices.  The  high  peak  powers  were  achieved  in  a series  of  experiments  us- 
ing intense  relativistic  electron  beams  (V  — I MV,  /~  30  kA)  and  powers  achieved  include 
900  MW  at  X - 4 cm'',  350  MW  at  X “2cm*,  8 MW  at  X - 8mm^and  2 MW  at 
X “ 4 mm®; however,  efficiency  of  converting  electron  beam  energy  to  EM  radiation  was  only 
~ 1%.  On  the  other  hand,  the  high  average  powers  were  generated  with  high  efficiency  and 
with  current  and  voltage  levels  similar  to  those  in  conventional  microwave  tubes;  thus,  the 


high  average  power  work  leads  directly  to  practical  devices. 


The  initiative  in  developing  high  average  power  gyrotron  was  first  taken  by  workers  at 
the  Gorkii  State  University  (USSR)'®".  The  key  to  achieving  efficient  devices  was  careful 
design  of  the  electron  gun.  In  the  Gorkii  studies,  a crossed  field  "magnetron  injection  gun"  was 
used  to  launch  an  annular  electron  beam  with  large  fraction  of  energy  transverse  to  the  axis 
and  with  minimum  energy  spread.  Work  on  nonuniform  cross-section  open  resonators  to  op- 
timize beam-wave  coupling  has  also  taken  place".  All  together  these  developments  have  led 
to  demonstration  of  a technological  breakthrough  in  CW  millimeter  wave  generation.  Using  a 
superconducting  magnet,  Zaytsev,  et.  a/.'®  have  generated  the  following  CW  powers;  12  kW  at 
X — 2.78  mm  with  31%  efficiency;  2.4  kW  at  X — 1.91  mm  with  9.5%  efficiency;  and  1.5  kW  at 
X “ 0.92  mm  with  6.2%  efficiency.  Other  impressive  results  by  Soviet  workers  include  the 
operation  of  a second  harmonic  device  in  a magnetic  field  of  only  6 kG  which  produced  at  X «• 
9 mm  a CW  power  of  10  kW  with  40%  efficiency  ".  Figure  1 compares  these  results  with  CW 
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power  available  from  other  microwave  devices;  an  advance  in  CW  power  capability  by  3 to  5 
orders  of  magnitude  is  clearly  indicated. 

The  Soviet  gyrotrons  have  operated  as  oscillators  using  a single  resonant  cavity  as  the  r.f. 
structure.  However,  the  gyrotron  process  lends  itself  to  the  use  of  a variety  of  r.f.  structures*, 
and  can  operate  as  an  amplifier  as  well  as  an  oscillator.  In  Fig.  2 schematics  of  three  types  of 
gyrotrons  are  sketched.  The  gyromonotron  in  Fig.  2a  corresponds  to  the  Soviet  oscillator 
configuration.  Figure  2b  shows  the  2-cavity  gyroklystron  amplifier  which  induces  a transverse 
phase  bunching  in  an  input  cavity,  allows  the  bunching  to  continue  ballistically  between  cavi- 
ties and  then  converts  transverse  electron  energy  to  wave  energy  in  an  output  cavity.  Opera- 
tion of  such  a device  was  reported  in  an  early  experiment  by  Wachtrel  and  Hirshfield and 
recently  by  Jory 

A third  type  of  gyrotron,  on  which  this  paper  focuses,  is  the  gyrotron  travelling  wave 
amplifier  (gyro-TWA,  Fig.  2c).  This  device  is  especially  uncomplicated  in  that  the  r.f.  structure 
is  a simple  waveguide.  It  has  a clear  advantage  in  handling  high  power  since  the  field 
strengths  encountered  for  a given  power  level  will  be  much  lower  in  a travelling  wave  device 
than  in  a device  employing  resonant  cavaties.  Also,  superior  bandwidth  and  tunability  charac- 
teristics are  expected.  The  operation  of  a gyro-travelling  wave  amplifier  was  demonstrated  in 
an  experiment  using  an  intense  relativistic  electron  beam'^.  A gain  of  16  dB  (1.1  dB/cm)  was 
achieved  at  8.6  GHz  with  a bandwidth  of  ~ 5%  and  an  output  power  of  4 MW;  there  was  also 
indication  of  octave-like  tunability  through  changing  the  magnetic  field  in  step  with  the  input 
frequency.  However,  efficiency  in  this  first  gyro-TWA  was  low  in  part  because  of  poor  beam 
quality. 

The  present  paper  describes  the  design  of  a millimeter-wave  gyro-TWA  which  is  intend- 
ed to  be  highly  efficient.  It  is  proposed  to  couple  the  r.f.  structure  designed  here  with  a "mag- 
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netron  injection  gun"  of  the  type  used  in  the  Gorkii  studies  so  than  an  electron  beam  can  be 
provided  with  sufficiently  high  quality  to  satisfy  the  design  assumptions. 


II.  DESIGN  CRITERIA  AND  THEORETICAL  APPROACHES 


The  methods  presented  here  for  parameter  optimization  are  generally  applicable  to  gyro- 
tron  travelling  wave  amplifiers.  As  an  illustration  we  will  consider  a particular  example,  name- 
ly, a gyrotron  travelling  wave  amplifier  aimed  at  the  laboratory  production  of  340  kW  (CW  or 
pulsed)  power  at  a frequency  of  35  GHz.  Wave  amplification  results  from  the  interaction 
between  the  circular  waveguide  mode. 


a»  ^ ~ O)}  “0, 

and  the  fundamental  cyclotron  harmonic  of  the  beam. 


(1) 


I “ ~ l^z^z  ~ " 0- 
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j where  cu,  is  the  cutoff  frequency  of  the  mode,  = ilf/y,  fl,  = cB/mc,  and  y = 

1 

-1-  v/)/c^l 
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i The  principal  design  criterion  is  to  attain  high  efficiency.  To  achieve  this,  it  is  desirable 

to  select  the  magnetic  field  such  that  the  two  curves  represented  by  Eqs.  (1)  and  (2)  intersect 
at  a grazing  or  near  grazing  angle  in  the  (a>,  k^)  plane,  i.e.  with  the  group  velocity  (k^c^/w)  of 
the  waveguide  mode  nearly  equal  to  the  beam  velocity  (v^).  The  reason  for  such  a choice  will 
become  evident. 

Figure  3 shows  the  configuration  of  the  beam-waveguide  system  on  which  our  calcula- 
tions are  based.  It  consists  of  an  annular  electron  beam  propagating  with  uniform  velocity  v, 
in  a circular  cross  section  waveguide  of  radius  r^.  the  system  is  immersed  in  an  external  mag- 
netic field,  Bq  ej,  oriented  along  the  axis  of  the  waveguide.  We  assume  that  all  electrons  have 
the  same  perpendicular  velocity  v^^.  Under  the  influence  of  the  external  magnetic  field,  the 


electrons  gyrate  at  the  frequency  (1^  with  the  Larmor  radius  We  assume  that 


the  guiding  centers  of  all  electrons  are  uniformly  distributed  on  a surface  of  constant  radius  rQ. 
Thus  the  beam  inner  radius  (r, ) and  outer  radius  (^2)  are  Ai  = /q  - and  = rb  + ri^. 

In  this  model,  we  have  assumed  a thin  (thickness  = 2 ri)  monoenergetic  electron  beam. 
Our  recent  computer  simulation  of  a carefully  designed  electron  gun  shows  that  a beam  with 
thickness  < 2.3  and  energy  spread  < < 1%  can  be  achieved  (see  Sec.  V). 


Calculations  of  the  optimum  operating  parameters  have  been  based  on  two  independent 
approaches.  First,  a linear  analytical  treatment  of  the  beam-wave  interaction  has  been  used  to 
calculate  the  linear  wave  frequency,  growth  rate,  and  the  beam  position  and  wave  number  for 
maximum  beam-wave  coupling.  As  in  earlier  analyses'^  the  presence  of  the  beam  is  treat- 
ed as  a perturbation  to  the  waveguide.  However,  the  problem  has  now  been  formulated  in 
realistic  cylindrical  geometry  rather  than  the  idealized  parallel  plate  geometry  previously  adopt- 
ed. This  improvement  in  modelling  is  of  considerable  importance,  since  efficiency  is  a sensi- 
tive function  of  quantities  which  depend  on  geometrical  factors.  A second  improvement  in 
modelling  is  that  we  treat  the  electron  guiding  center  radius  zq  as  a variable  parameter,  while  in 
earlier  parallel-plate  models'^  the  guiding  centers  were  centrally  located  between  the 
plates.  This  gives  us  the  means  to  optimize  the  beam-wave  coupling  with  respect  to  rg. 


For  generality,  all  the  calculations  have  been  done  in  the  beam  frame  (i.e.  the  reference 
frame  in  which  the  beam  axial  velocity  is  zero).  Furthermore,  the  following  nomalization  pro- 
cedures are  introduced,  through  which  the  waveguide  radius  is  scaled  out  of  the  problem: 


length  normalized  to  r„ 
frequency  normalized  to  dr^. 
velocity  normalized  to  c 


(3) 

(4) 

(5) 


5 


I 

f 

i 


i 


i 


EM  field  normalized  to  mc^ ler„. 


P 

L 


(6) 


The  results  obtained  are  thus  applicable  to  any  waveguide  radius  and  for  any  beam  ax- 
ial velocity  v^  by  the  appropriate  beam  to  lab  frame  transformation.  In  Appendix  A we  have 
included  the  transformation  formulae  for  conversion  from  beam  frame  to  lab  frame  quantities 
and  from  dimensionless  to  physical  units. 


To  maintain  a clear  distinction  between  lab  frame,  beam  frame,  normalized  and  physical 
quantities  the  following  notation  has  been  adopted: 


(i)  Beam  frame  quantities  are  denoted  by  primed  symbols  (e.g.,  -y'),  laboratory  frame 
quantities  are  unprimed  (e.g.  y). 


(ii)  Normalized  dimensionless  quantities  are  denoted  by  a bar  (e.g.  flj,).  Unbarred  sym- 
bols represent  either  physical  quantities  (e.g.  flj)  or  naturally  dimensionless  ones  (e.g.  y). 


The  derivation  of  the  linear  dispersion  relation  is  lengthy  but  straightward.  Since  the 
methods  are  analogous  to  those  described  in  detail  by  Ott  and  Manheimer, only  the  result 
will  be  presented  here. 


The  normalized  dispersion  relation  for  the  excitation  of  the  mode  at  the  fundamen- 
tal electron  cyclotron  harmonic  (expressed  in  beam  frame)  is 


w'2 


4i/' 


y'J$  (x,  ) 


(co  H (x,  Fq  ) 


(w' 


w'(?(x,?o,jf,ri) 


r - n. 


(7) 


where  Xj  —3.832  is  the  first  nonzero  root  of  7]  (x)  — 0,v'=N'rf  is  a dimensionless  beam  den- 
sity parameter  (N' is  the  total  number  of  electrons  per  unit  length  r,=eVmc^  is  the  classical 


electron  radius),  = Vj^/c,  and 
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^ . H(x^)  = [j,  (x)/)  (>)]^ 

t 

[ Q(.x^)  = 2Hix^)  +>'7,  (>’)7i  (^)U(^  (x)  (1  +x  ~^)  + [7,  (x)l^) 

[ t 

j + 27,  (x)7,  (x)7,'  (>>)  [>7,  (;»)  - 7,  (>')  Mxy 


(8) 


(9) 


Note  that  the  primes  in  Eqs  (8)  and  (9)  represent,  respectively,  the  first  and  second  derivatives 
of  the  Bessel  function  7,  with  respect  to  its  argument.  Elsewhere  else  in  this  paper,  the  primes 
indicate  beam  frame  quantity. 


In  Eq.  (7),  and  are  frame  independent  quantities.  It  can  be  easily  shown  that 
the  normalized  cutoff  frequency  of  the  TEq,  mode,  aT^.,  is  numerically  equal  to  x, . 

This  analytical  work  has  been  complemented  by  a single  wave  numerical  simulation  code. 
The  Cartesian  version  of  this  code  was  initially  developed  by  Sprangle  and  Manheimer. It 
was  later  refined  by  Sprangle  and  Drobot’^  and  employed  to  investigate  the  saturation  mechan- 
isms. A complete  documentation  of  the  Cartesian  code  appears  in  Ref.  17.  Here  we  have  con- 
verted the  code  into  cylindrical  geometry  and  use  it  to  simulate  the  beam-wave  interaction 
processes  until  the  wave  reaches  saturation.  A self  consistent  value  of  the  efllciencv  can  there- 
fore be  determined.  The  simulation  shows  two  different  but  simultaneously  present  saturation 
mechanisms  —depletion  of  free  electron  energy  and  electron  trapping  in  phase  space.  For  a de- 
tailed description  of  these  mechanisms,  see  Ref.  17.  The  values  quoted  for  the  efficiency  in 
Sec.  Ill  are  those  taken  from  the  simulation  data.  Conservation  of  energy  was  monitored  in  the 
code  and  never  deviated  by  more  than  0.1%.  As  a cross  check,  the  linear  frequency  and  growth 
rate  obtained  from  Eq.  (7)  have  been  compared  to  the  linear  stages  of  the  simulation.  The 
agreement  between  the  two  approaches  was  very  good,  (see  Sec.  III). 
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III.  BEAM  FRAME  PARAMETER  OPTIMIZATION 


I 


As  seen  from  Eq.  (7),  solutions  of  the  problem  in  the  beam  frame  and  in  dimensionless 
units  require  the  specification  vc  parameters:  ,y  , and  11^.  In  the  present  design, 

we  will  choose  and  k.  so  that  the  beam-wave  coupling,  represented  by  the  linear  growth  rate 
w, , is  a maximum  with  respect  to  these  two  parameters.  The  parameter  </'wiil  be  determined 
from  the  required  output  power.  With  these  three  parameters  specified,  the  energy  conversion 
efficiency  i)',  defined  as  the  ratio  of  the  final  wave  energy  to  the  initial  beam  energy,  can  be 
calculated  numerically  varying  the  two  remaining  parameters  -y' and  II  j.  The  point  in  the 
{y  ,llj ) parameter  space  where  tj' peaks  will  then  be  chosen  for  the  design. 

The  first  term  on  the  RHS  of  Eq.  (7),  proportional  to  is  the  driving  term  for  the  ins- 
tability. Hence,  to  maximize  the  beam-wave  coupling  with  respect  to  Tq,  we  let  Tq  - 0.48  so 
that  peaks.  This  is  also  the  radius  at  which  the  wave  electric  field 

Eu  (a  J\  Ujr)]  peaks.  We  mention  in  passing  that  for  operation  at  higher  cyclotron  harmon- 
ics, the  optimum  beam  position  does  not  in  general  coincide  with  the  peak  of  the  wave  electric 
field. 

The  second  term  on  the  RHS  of  Eq.  (7)  imposes  a threshold  energy  on  the  beam,  below 
which  there  will  be  no  instability.  For  most  gyrotron  devices  including  the  presently  designed 
one,  the  beam  is  well  above  threshold  energy.  Hence  this  term  can  be  neglected.  We  expand 
ilfj  about  a fixed  value  expand  w'and  k^  about  the  intersecting  point  (tog.  ^ 

the  two  equations. 


Substituting 


CO  “ tOQ  + Aco 


“ ^bO  + 

into  Eq.  (7)  and  keep  terms  to  first  order  in  AoT  , A/tj',  and  AfT*,  we  obtain 


A-'  '‘^0  *r' 

a<ti  — —7—  A/c^ 

coq 


(Aco  -All*)  = . , . 

y Jq  (jr,  )coo 


(12) 


In  Eq.  (12),  Aco  >=  Aco^  + /Aa>,  is  a complex  number,  while  all  other  quantities  are  real. 
Differentiating  Eq.  (12)  with  respect  to  AiL,  we  obtain 


d{Lto  ) ^ ~ ^ ^ i ) 

d{^k^)  (oq  (3Aco  ~A(1j  ~ 2k,Q Akj/tuQ ) (13) 

Eq.  (13)  shows  that  if  /c^q  «*  0,  the  growth  rate  Aco,  is  a maximum  regardless  of  the 
values  of  Ak^  and  AD*,  i.e.,  Aco,  remains  a maximum  with  respect  to  k,  even  with  small  varia- 
tions in  n*.  In  comparison,  if  k^^  ^ 0,Aw,  peaks  [i.e.  the  RHS  of  Eq.  (13)  is  real]  only  when 
^^2  “^(^6  “ 0-  The  case  /c,q  =0  thus  has  the  advantage  of  insuring  maximum  growth 
rate  (with  respect  to  F^)  when  the  magnetic  field  is  varied  to  optimize  the  efficiency.  It  implies 
a grazing  intersection  of  the  two  curves  represented  by  Eqs.  (10)  and  (11)  (see  Fig.  4).  This  ex- 
plains why  we  have  chosen  a grazing  or  near  grazing  intersection  for  the  design. 


It  is  convenient  to  define  a new  parameter  X'\ 

X = n*7co^. 

When  the  applied  magnetic  field  is  such  that  A"=  1 the  two  curves  (Eq.  (10)  and  (11)]  intersect  at  an 
exact  grazing  angle. 
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We  now  consider  the  choice  of  v For  each  value  of  v one  can  calculate  the  point  in  (y 
space  where  r?'  peaks.  Given  the  beam  axial  velocity  v^  in  the  lab  frame,  the  beam  power  (Z*^)  and  the 
output  wave  power  at  the  point  of  maximum  efficiency  can  then  be  calculated  (see  Appendix  A, 
Eqs.  (A5-7)] . To  specify  v^,  we  define  a lab  frame  parameter  a, 

a = yjv,. 

Thus,  for  a fixed  value  of  a,  there  is  a one  to  one  correspondence  between  and  v . 
Lengthy  calculations  are  required  to  deduce  p from  P^.  because  one  might  have  to  scan  a 
range  of  values  for  p before  finding  the  one  that  yields  the  prescribed  P^.  In  practice,  howev- 
er, one  can  make  a close  initial  bucss  for  p ' on  the  basis  of  the  required  output  wave  power, 
the  expected  efficiency  and  beam  energy  [see  Eqs.  (A5-7)l,  and  then  proceed  to  calculate  the 
optimum  y ',  Jf , and  P^.  If  P^  turns  out  to  be  too  large,  for  example,  by  10%,  then  the  next 
(iuess  for  v \ 10%  smaller,  can  be  very  accurate  because  the  relation  between  P^  and  p is  al- 
most linear.  If  there  is  an  acceptable  range  of  P^,  the  initial  guess  of  p usually  suffices. 

The  choice  for  the  velocity  ratio  a is  guided  by  two  considerations  - having  as  much  en- 
ergy in  the  perpendicular  gyromotion  as  possible  and  minimizing  the  beam  temperature  spread, 
ihe  second  consideration  sets  an  upper  limit  on  a.  Our  electron  gun  simulation  code  showed 
that  for  a > 1.5,  the  beam  temperature  began  to  increase  rapidly  and  was  no  longer  suitable. 
Thus  we  let  « — 1.5.  The  required  output  wave  power  is  ~ 300  kW,  which  corresponds  to  p 
= 0.002,  if  we  assume  y =1.1  and  a 45%  efficiency.  Thus  we  let  v “0.002  be  our  initial 
guess.  We  now  present  the  calculations  that  lead  to  the  optimized  y , X , and  for  the  ini- 
tially chosen  p '.  The  calculated  optimum  P^  (340  kW)  turns  out  to  be  about  11%  larger  than 
expected,  but  falls  in  the  acceptable  range. 

With  p'  specified  at  p'  - 0.002  the  efficiency  tj'  has  been  calculated  in  a two  parameter 
space  (y  ,X  ).  The  results  are  summarized  in  Figures  5 and  6.  Fig.  5 shows  two  plots  of  tj'  vs 
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y'.  The  solid  curve  represents  the  efficiency  for  grazing  interaction  (A"  =1).  It  predicts  a 
maximum  efficiency  tj'  = 42%  at  y’  — 1.04.  The  dashed  curve  shows  the  efficiency  as  a func- 
tion of  y but  for  an  optimized  value  of  X'.  The  peak  efficiency  is  now  tj'  = 73%  at  y'  = 1.10 
and  X'  = 0.9S7.  It  is  interesting  to  note  that  the  optimization  in  X'  shifts  the  point  of  max- 
imum efficiency  to  higher  y'  than  found  for  grazing  interaction.  Decreasing  X'  from  1 has  the 
effect  of  delaying  phase  trapping'^,  which  is  one  of  the  two  mechanisms  responsible  for  satura- 
tion. Since  phase  trapping  is  more  dominant  at  higher  y'  this  qualitatively  explains  why  the 
efficiency  peak  moves  to  higher  y'  with  magnetic  mismatching.  To  determine  quantatively  the 
value  of  y'  where  the  efficiency  peaks,  a thorough  search  of  the  (y'  , X')  parameter  space  was 
required. 

It  can  be  seen  in  Fig.  S that  at  the  value  of  y'  for  which  the  maximum  efficiency  was 
found  a small  variation  in  if'from  1.0  to  0.96  caused  a large  change  in  the  efficiency  from  30% 
to  70%.  To  further  illustrate  the  sensitivity  of  the  efficiency  as  a function  of  A"  a graph  of  17 'vs 
X'  is  shown  in  Fig.  6 for  y'  — 1.10,  the  value  corresponding  to  peak  efficiency.  It  can  be  seen 
in  Fig.  6 that  tj'  is  a sensitive  function  of  X'.  It  increases  sharply  as  A”  is  decreased  from  unity 
and  then  falls  precipitously  if  AT' is  reduced  too  much.  This  occurs  because  for  too  low  a value 
of  AT',  the  system  becomes  linearly  stable.  We  would  like  to  stress  that  for  an  experiment  to 
operate  near  the  peak  value  of  t)',  it  is  essential  to  control  parameters  such  as  magnetic  held, 
beam  voltage  and  driver  frequency  within  an  accuracy  of  <1%. 

Fig.  7 shows  the  growth  rate  (ci>,  ) and  the  difference  between  wave  frequency  ) and 
electron  cyclotron  frequency  (H^)  as  a function  of  A”  for  the  same  parameters  is  used  in  Fig. 
6.  Solid  curves  are  calculated  from  Eq.  (7)  and  the  solid  data  points  are  taken  fiom  the  simula- 
tion runs.  The  agreement  between  analytical  and  simulation  results  is  found  to  be  very  good 
for  the  runs  shown  in  this  figure  as  well  as  for  all  other  runs.  From  Figs.  6 and  7,  we  observe 
that  when  the  efficiency  approaches  a maximum  as  a result  of  varying  X’,  the  growth  rate  ap- 
proaches zero.  Too  low  a growth  rate  is  clearly  undesirable.  However,  as  one  moves  slightly 
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away  from  the  peak  efficiency,  the  growth  rate  increases  rapidly.  Some  compromise  must  be 
made  between  growth  rate  and  efficiency  in  selecting  the  operating  point.  For  the  designed 
amplifier,  we  have  chosen  the  point  X'  * 0.96  where  the  efficiency  is  n “ 70%,  only  three  per- 
cent below  the  peak  value.  As  will  be  shown,  this  point  corresponds  to  a reasonably  high 
power  gain  (2  dB/cm).  The  other  parameters  corresponding  to  this  point  arc:  v -0.002,7^  - 
0.48,  y ' - 1.1,  k’j  -0,  w,'  - 3.791,  and  w,'  -0.0342. 

IV.  CONVERSION  OF  BEAM  FRAME  OPTIMUM  PARAMETERS 

INTO  DESIGN  PARAMETERS 

To  convert  the  beam  frame  optimum  parameters  into  actual  design  parameters,  two  peices 
of  information  are  required: 

(i)  the  beam  axial  velocity  (in  lab  frame),  which  can  be  calculated  from  the  previously 
specified  value  of  a. 

(ii)  the  waveguide  radius,  which  is  determined  from  the  desired  wave  frequency,  / - 35 
GHz.  [see  Appendix  A,  Eq.  (A-4)I. 

Using  the  conversion  formulae  tabulated  in  Appendix  A,  we  have  converted  the  op- 
timum parameters  obtained  in  Sec.  II  into  the  design  parameters  for  a 35  GHz,  340  kW  gyro- 
TWA.  They  are  shown  in  Table  I. 

The  experimental  configuration  which  was  based  on  these  design  parameters  is  shown  in 
Fig.  8.  The  entire  device  is  placed  within  the  bore  of  a superconducting  magnet.  This  allows 
precise  shaping  of  the  magnetic  field  over  the  entire  electron  beam.  Such  control  is  essential 
because  of  the  great  impact  small  changes  in  the  magnetic  field  have  on  both  the  device 
efficiency  and  the  operation  of  the  electron  gun.  The  magnetic  field  profile  shows  a carefully 
designed  compression  region  required  to  obtain  a beam  with  the  proper  velocity  ratio  a.  After 


12 


being  compressed,  the  beam  enters  a uniform  magnetic  field  region  where  the  cyclotron  maser 
mechanism  operates.  An  external  driver  generates  a r.f.  signal  which  is  injected  at  two  azimu- 
thal positions  at  the  left  end  of  the  drift  tube  so  as  to  launch  a TEgi  wave.  The  use  of  an  ab- 
sorber ensures  that  the  wave  can  propagate  in  only  one  direction  so  as  to  prevent  oscillations. 
Amplification  oc  .urs  in  the  uniform  field  region  to  the  right  of  the  r.f.  injection  point.  Beyond 
this  region  the  electron  beam  is  collected  at  the  wall  as  the  magnetic  field  falls  off.  Stray  elec- 
trons are  deflected  by  a small  magnet.  The  rJ.  signal  continues  down  the  drift  tube  and  exits 
at  the  window. 

V.  DISCUSSION  AND  CONCLUSION 


At  this  time  there  exists  no  comprehensive  theory  detailing  the  effects  of  b'.am  tempera- 
ture on  the  efficiency  of  gyro- TWA ’s.  It  is  possible,  however,  on  the  basis  of  nonlinear  satura- 
tion mechanisms  to  derive  the  following  qualitative  condition  for  the  validity  of  the  cold  beam 
approximation 


Sy  '/y  <<  2(y  - I )7)/y, 

(14) 

where  fi-y  is  the  energy  spread  in  the  beam  frame  and  t)  is  the  laboratory  efficiency  calculated 
for  a cold  beam.  A numerical  simulation  has  been  made  of  the  NRL  gun  (see  Fig.  8),  which  is 
similar  in  configuration  to  the  electron  gun  used  by  Kisel**  et  al.  The  beam  produced  in  the 
simulation  is  almost  monoenergetic  (typically,  Sy/  iy  — 1 ) < < 1 %1  so  that  the  energy 

spread  in  the  beam  frame  (8y  ) is  caused  mainly  by  a velocity  pitch  angle  spread  S0 
(0  =tan  a)  of  the  beam  electrons.  For  the  beam  parameters  in  Table  I,  we  find  that  80  = 
0.049  and  6y/  {y  — 1 ) = 0.002.  When  the  beam  energy  spread  (in  lab  frame)  is  much  small- 
er than  pitch  angle  spread,  Eq.  (14)  can  be  replaced  by 


80  < < 


2 (y  — 1 )'r}  (1  - p}  ) 

yJTK 


(15) 


where  = yjc  and  = v^/c. 
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Substituting  the  beam  parameters  in  Table  I into  Eq.  (1 S)  gives 

< < 2.1  T»  . 

This  is  well  satisfied  since  60  ^ 0.0S  and  the  projected  value  for  t)  is  — O.S.  The  cold  beam  ap- 
proximation adopted  in  our  theory  is  therefore  realistic  based  on  the  simulation  studies  of  the 
electron  gun. 

In  conclusion,  we  have  formulated  the  linear  cyclotron  maser  theory  and  a numerical 
simulation  code  in  cylindrical  geometry  and  shown,  through  an  example,  how  all  the  parame- 
ters of  a gyrotron  device  are  optimized  to  achieve  high  efficiencies.  The  choice  of  beam  posi- 
tion yg  and  wave  number  are  based  on  the  condition  of  maximum  linear  beam-wave  cou- 
pling. The  beam  density  parameter  v is  determined  from  the  required  power  level.  Finally, 
the  specification  of  the  lab  frame  beam  energy  y and  the  cyclotron  frequency  parameter  is 
based  on  extensive  numerical  calculation  of  efficiency  in  the  (y  ,X  ) parameter  space. 

Sensitivity  of  the  efficiency  to  small  variations  in  the  relativistic  electron  cyclotron  fre- 
quency requires  that  the  beam  energy  and  the  external  magnetic  field  be  controllable  within  an 
accuracy  of  < 1%. 

The  authors  would  like  to  thank  Dr.  T.  Godlove,  Dr.  R.  Parker,  Dr.  M.  Reed,  and  Dr.  P. 


Sprangle  for  many  helpful  discussions. 
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Appendix  A 

CONVERSION  FORMULAE 


Theoretical  calculations  in  the  present  work  are  carried  out  in  beam  frame  and  parame- 
ters are  normalized  according  to  Eqs.  (3)-(6).  In  this  appendix,  we  tabulate  the  formulae  need- 
ed to  convert  beam  frame  normalized  quantities  to  lab  frame  phyical  quantities.  Column  2 of 
Table  II  converts  the  beam  frame  normalized  quantities  in  column  1 to  lab  frame  normalized 
quantities.  The  information  needed  for  the  conversion  is  the  axial  velocity  v^  (in  lab  frame) 
which  defines  the  quantity  in  column  2, 

y,=  [1  —\}/c^] 


Formulae  presented  in  column  2 are  based  on  the  Lorentz  transform.  Derivations  are  ob- 
vious except  for  the  following  items. 

(/)  T)  - (y' -1  )/(y^y' — 1 ) 

(A.l) 

To  derive  Eq.  (A.l),  we  note  that  the  time  averaged  EM  field  energy  per  unit  length  is 
(for  TE„„  mode) 


- 


where  A is  a geometrical  factor,  is  the  axial  wave  magnetic  field  amplitude,  and  is  the 
waveguide  cutoff  frequency. 


The  injected  beam  energy  per  unit  length  is 


"N(y-l)mc^ 


where  S is  defined  in  the  text. 


The  lab  frame  efficiency  (t) ) can  therefore  be  written 


IS 


<Wf> 


K B}  me  W 


ina}jN(y  — 1 ) 

In  Eq.  (A.2)  K,  Bj  and  otg  are  all  frame  independent  quantities. 


(A.2) 


The  beam  frame  efficiency  (t)  ')  is  also  expressed  by  Eq.  (A.2)  with  u,  N,  and  y replaced 
by  at',  N',  and  y ' respectively.  Thus 


jL  a,^\'(y'  -1) 
V'  w'^N(y-l) 


(A.3) 


Since  <u  - y^w',  y - and  N - y^N',  substitution  of  these  relations  into  Eq.  (A.3) 
yields  Eq.  (A.l). 


(//)  (Of  “ y^w^'and  w,-  “(!>,■ '/y^ 

In  the  beam  frame,  —0.  Hence  all  quantities  vary  as  exp  ( —i<o'i').  To  transform  this 
expression  into  the  lab  frame,  we  use  the  relation  t’  — y^  (t  — v^z/c^^  and  obtain 

exp  ( —i(o '/')  “ exp  ( — / (w, ' + i<o ')  y ^ (/  — \fZ/c^ ) ] 


-exp 




-/y^O),  t + I 


- . yz^i'^z 

j 2 + y^toj  t — z 

cZ  f.1 


Thus,  the  exp(  —i<o'i')  dependence,  viewed  from  the  lab  frame,  has  the  following  wave 
frequency  (w,.)  and  wave  number  (k^ ), 


kf  -y^(0,V^/c^ 

and,  if  one  follows  a beam  segment  (i.e.  let  z — v^r),  one  observes  a growth  rate  given  by 

‘"/  - y^a>,  '(l  -v/Zc^)  -w,  '/y^. 

Column  3 converts  the  lab  frame  normalized  quantities  in  column  2 to  lab  frame  physical 


quantities.  The  conversion  formulae  are  based  on  Eqs.  (3)- (5)  in  the  main  text.  The  informa- 
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tion  needed  at  this  step  is  the  waveguide  radius  which  can  be  determined  from  the  desired 
wave  frequency  through 


(A.4) 

where  uT,  is  the  calculated  wave  frequency  (normalized)  and  / is  the  desired  wave  frequency  in 
Hz. 

From  the  lab  frame  physical  quantities  in  column  3,  we  can  derive  the  beam  current 
(/{,),  beam  power  (Pj,),  the  maximum  output  wave  power  (P^),  and  the  power  gain  G : 


/j  ^pxnc^vje  — 1.707  X 10^  vv^/c  Amp, 

Pj  ”■ /j,  (y  — l)mc^  — 8.535  X 10®  V (y  — l)v^/c  kW, 

- n Pb> 

G — 8.7w,/Vj  dB/unit  length. 


(A.5) 

(A.6) 

(A.7) 

(A.8) 
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Table  I — Efficiency  optimized  operating  parameters 
for  the  NRL  gyrotron  travelling  wave  amplifier, 
a = 1.5,  / = 35  GHz  (X  = 8.57  mm). 


V (density  parameter) 

2.076  X 10-3 

V (electron  energy) 

70.82  keV 

If,  (beam  current) 

9.48  Amp 

r}  (effieciency) 

51.0  % 

Pf,  (beam  power) 

671.5  kW 

(wave  power) 

342.5  kW 

Bo 

12.87  kG 

K 

1.96  cm'l 

5.37  mm 

'•o 

2.52  mm 

0.61  mm 

'•l 

1.91  mm 

'‘2 

3.13  mm 

vi/c 

0.401 

vjc 

0.268 

G (power  gain) 

2.0  dB/cm 

Table  II  — Formulae  for  transformation  from  beam  frame  into  lab  frame  and 
for  converting  normalized  quantities  into  physical  quantities 


Beam  Frame 
Normalized  Quantities 

Lab  Frame 

Normalized  Quantities 

Lab  Frame 
Physical  Quantities 

p' 

11 

p (dimensionless) 

y 

7 = 72  7' 

7 (dimensionless) 

n' 

V = t?'72(7'-1)/(727'-1) 

17  (dimensionless) 

“r  = 72  w' 

(Or  = COrCAu, 

CO- 

w,-  = ^■/7^ 

iOf  = (0,c/r^ 

*2  = 72  tJr 

K = k^/r^ 

n,=n' 

Slg  = cAu, 

r'L 

ri  = 

Tl  = ryr^ 

^2 

^2  = ^2 

’■2  = 

20 


Fig.  3 - CTOss^ectional  view  of  the  beam  modeL  The  applied 
magnetic  field  (not  shown)  points  toward  the  reader.  The 
electrons  are  monoenergetic  and  all  have  the  same  Larmor 
radius  r^.  Guiding  centers  of  all  electrons  are  uniformly 
distributed  on  the  circle  of  constant  radius  r^. 


Fig.  4 - The  vacuum  wave  guide  characteristic  curve  (Eq. 
£10)1  and  the  beam  characteristic  curve  [Eq.  (11)  in  the  (w', 


k,)  plane. 
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Fig.  5 - t)'  venvu  7'  for  grazing  interMction  case  (.X'  = 
1,  solid  curve)  and  for  optimized  near  grazing  inter- 
section case  (X'  < 1,  dashed  curve).  A point  on  the 
dashed  curve,  tj'  (y',X'),  is  obtained  by  varying  X'  but 
keeping  y fixed  until  tj'  reaches  a maximum, 

illustrated  in  Fig.  6.  The  peak  of  the  dashed  curve  is 
therefore  the  peak  of  tj'  in  the  (y',X)  parameter  space. 
The  parameters  used  for  this  figure  are  = 0.48,  i’  = 
0.002,  and  if  ^ = 0. 


Fig.  6 - Efficiency  tj  versus  Jf 'for  y • i.l,  0.48,  t>' 
0.002,  and  - 0. 


Fig.  8 - Conflguiation  of  the  designed  gyro  travelling  wave  amplifier. 


